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I. INTRODUCTION 



Because of dimensional transmutation, the properties of SU(N) gauge theories are of 
great interest. Once one mass scale is set, all other properties are in principle determined. 
For example, if the theory exhibits a deconfining phase transition at a temperature T^, then 
the nature of the phase transition is completely determined once the value of is known. 

Understanding the deconfining phase transition at small N requires numerical simula- 
tions on the lattice [1, 2]. Using these results, a matrix model was developed to model 
deconfinement. These models involve zero [3, 4], one [5], and two [6, 7] parameters. They 
are soluble analytically for two and three colors, and numerically for four or more colors. 

This matrix model is also soluble in the limit of infinite N [8]. The phase transition at 
infinite N is exceptional, and can be termed a Gross- Witten-Wadia transition [9-11]. For the 
deconfining transition, such a transition was first seen to occur on a femto-sphere [12-16]. 
For such a transition, at infinite N it exhibits aspects of both first and second order phase 
transitions. It is of first order in that the latent heat is nonzero and proportional to ~ N 2 . 
It is also first order in that the Polyakov loop jumps from to \ at On the other hand, 
it is of second order in that several quantities, such as the specific heat, exhibit nontrivial 
critical exponents. Such an unusual transition only occurs at infinite N, as at finite N the 
transition is of first order. 

Consider a phase transition which can be characterized by the change of a single field. If 
the transition is of first order, then at the transition temperature there are two degenerate 
minima, with a nonzero barrier between them. Then the order-disorder interface tension 
is nonzero, given by the probability to tunnel between the two minima. Conversely, if the 
transition is of second order, at the transition temperature the two minima coincide. Then 
there is no barrier to tunnel between them, and the order- disorder interface tension vanishes. 

In a gauge theory asking about the potential is more subtle. In a matrix model of an 
SU(N) gauge theory, the potential exists in N — 1 dimensions (the number of commuting 
diagonal generators). A physical question is to ask how the interface tensions behave, es- 
pecially at the phase transition. There are two such interface tensions. There is the usual 
order-disorder interface tension, which exists only at T^. There are also order-order interface 
tensions [17-22] , which are directly related to the 't Hooft loops for Z(N) charges [23, 24]. 
These exist for all T > T^. 
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In this paper we consider the interface tensions in the matrix model near T^. We find that 
in the matrix model, both the order-order interface tension and the order-disorder interface 
tensions, vanish identically at T = Td when N = oo. 

There is a simple heuristic explanation for our results. Consider the potential for the 
simplest Polyakov loop, that in the fundamental representation, E± = (l/iV)trL, where L 
is the thermal Wilson line. Then the Gross- Witten-Wadia transition occurs because at 
infinite N, the potential for t\ is completely flat between and \. That is, at Td there 
are two distinict minima, as expected for a first order transition. Nevertheless, because the 
potential is flat at N = oo between the two minima, there is no barrier to tunnel between 
them, and so the interface tensions vanish. Such a flat potential was first found for the 
Gross- Witten-Wadia transition on a femto-sphere [12-16]. 

As noted, in a gauge theory there are other degrees of freedom. For example, one can 
consider higher powers higher powers of L, lj = jjtrV, for j = 2 . . . N - 1. For the Gross- 
Witten-Wadia transition, though, all £j vanish at Tj when j > 2, which allows us to consider 
the potential as a function of a single variable, t\. We stress, however, that this is true only 
at Td, and not for T ^ Td- Indeed, while we estimate the behavior of the (order-order) 
interface tensions for T > Td, we cannot solve the problem in full generality. This is because 
away from Td, all £j contribute, and it is a much more difficult problem than at T^. 

The order-disorder interface tension has been measured through numerical simulations 
on the lattice at iV = 3 by Beinlich, Karsch, and Peikert [25]. Lucini, Teper, and Wenger 
give results for the order-disorder interface tension for N = 3, 4, and 6, and use these to 
extrapolate to N = oo [26]. On general grounds the order-disorder interface tension should 
be proportional to ~ N 2 T%, with a coefficient which is naturally of order one. Instead, Ref. 
[26] finds a very small coefficient at iV = oo, ~ .014. In the matrix model the coefficient 
~ N 2 vanishes identically, and the true behavior is only ~ N. This is because in the matrix 
model, the height of the barrier between the two distinct minima is small, ~ 1/N. We 
discuss this further in Sec. (IV) and in the Conclusions, Sec. (V). 

II. REVIEW OF LARGE N THERMODYNAMICS 

We are interested in the thermodynamics of pure SU(N) gauge theory for temperatures 
a few times that for the deconfining phase transition at Td- The order parameter of the 
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thermodynamics is taken to be the Wilson line: 



L(x) = Pexp ^ig J A (x,T)dr^j . (1) 
By a gauge transformation, we can diagonalize the field A as 

AS = ^yn (2) 

where i, j = 1 ■ ■ ■ N and the eigenvalues g« are subject to the SU(N) constraint q t = 0. 
In our model are the fundamental variables to characterize the transition. We assume 
that after integrating out the other components of gluon field Ai, that we obtain an effective 
potential for qi [3-8]. 



V ef f(q) = + d 2 (T)V 2 (q), (3) 

N 

V n (q) = Y l \q i -q j \ n (l-\q i -q J \) n . (4) 

The potential includes both perturbative (V 2 ) and non-perturbative (Vi) contributions. The 
temperature dependent functions d\ and d 2 are given by 

di{T) = %c lT > Tl d 2 (T) = f (T 4 - c 2 T 2 Tj) . (5) 

At tree level, the kinetic term is 

In a mean-field approximation, the kinetic term does not contribution for a spatially homo- 
geneous states, and so it can be ignored for thermodynamic quantities. The kinetic term 
does enter in computing the interface tension in the following sections. 

In the infinite N limit we introduce a continuous variable x = j^. Labeling the eigenvalue 
qi —> q(x), we introduce the eigenvalue density 

N 



1 f 1 

p(q) = lim — V S(q - q { ) = I 

N^oo iV / — ' J 



dx5[q — q(x)] = — . (7) 



At finite N, the identities 

N N 
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become, at infinite N, 

J dqp{q) = 1, J dqp{q)q = 0. (9) 

The potential is proportional to iV 2 , 

V n (q) = N 2 V n (q) = N 2 J dx dy \q(x) - q(y)\ n (1 - \q(x) - q(y)\) n 

= N 2 Jdq dq' p(q) p(q') \q - q'\ n (1 - |g - q'\) n . (10) 

This representation transforms the potential into a polynomial in q. 
The minimum of Eq. (10) was found in Ref. [8]. The solution is 

p(q) = 1 + bcosdq, -q < q < q Q , (11a) 



W 1 ?- (nb) 

cot{dqo) = i{l- qo )-d(i/2-q y (llc) 

dWl V d 2 fl ^ 2 



-'/o) +yU~ 9o J +L (lld) 

For T > T d , d > 2n and q < h The eigenvalues do not span the full range between — | 
and |. The density is discontinuous at the end points p(±g ) > 0. For T = T£, qo = \ and 
d = 2ir, the density is continuous for all values of q in [—1/2, 1/2]. In particulalr it vanishes 
at the end points p(±go) = 0. For T < Td, the theory is in confined phase, with a uniform 
distribution of eigenvalues over the unit circle, Eq. (11a) with q = | and 6 = 0. 

For the potential Eq. (10) and the eigenvalue distribution Eq. (11a), at T d , q = |, and 
the potential is independent of b. Changing b from to 1 interpolates between confined and 
deconfined phase, but does not change the potential. Hence b is a zero mode of the potential 
corresponding to changing the overall shape of the distribution. This will play an important 
role in the construction of the interface. For short we call the change of shape related to b 
the fe-mode. 

It is also worth emphasizing that in the derivation of Ref. [8] , we have assumed that the 
eigenvalue density is symmetric in q. We can obtain different distributions by applying an 
arbitrary transform to a given solution. A transform of charge k, k = 1 . . . (N — 1), 
is given by 

k k k-N k-N , . 

<?i,<?2,- ■ -<?tv -» qi + —,---q N - k + j^,qN-k+i+ N ,---qw+ N ■ (12) 
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Assuming q± < q<i < ■ • • < q^ and \qi — q<j\ < 1, we can relabel the eigenvalues such that 
they are in an increasing order: 

k — N k — N k k 

qi,<h,--qN ^q N -k+i+ N ,-"Qn + N } qi + —,---q N -k + ^- (13) 

In the infinite N limit, the tranform takes the following form. Define the inverse 
function of x(q) = f_ o dq'p(q') as Q(x), 

|q(i+1-A)-1 + A 0<x<A 
q = Q( x ) q = 1 . (14) 

I Q(x - A) + A A<x<l 

Since the potential is invariant under Zn transformations, smooth changes in A are another 
zero mode of the potential, which we call the shift mode. This is also relevant to the 
construction of the interface. 

We stress that both the b and shift modes are become zero modes only at infinite N. The 
former is because of the flatness of the potential. The latter is because the Z^ symmetry 
becomes a continuous U(l) at N = oo. We comment on what happens at finite A^ later. 

III. INTERFACE TENSION 

An interface is a topological object interpolating between two vacua of the theory. Sup- 
pose the two vacua are separated in the z direction and extended in the x — y plane. The 
effective action has an area law when the transverse size L tr is large: S e ff = <xL tr . The 
proportionality constant defines the interface tension. Up to cubic order in the perturbative 
expansion, the order-order interface tension exhibits Casimir scaling [19, 20] 

aock(N-k). (15) 

An important question to address is whether the order-order interface tension satisfies 
Casimir scaling in the matrix model. 

An interface tension necessarily involves a spatial gradient along the z direction. Conse- 
quently, the kinetic term must be included. Since the potential is simple when written in 
terms of the eigenvalue density p(q), it is useful to write the kinetic term in terms of the 
same variable as well. Assuming that there is a spatial gradient only along the z direction, 

* (9) = m = N>Vfl j dzdx (M^!) 2 , (16) 
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where the partial derivative is taken at fixed x. 

We start with the eigenvalue density p(q, z) = dx(q, z)/dq. Assuming that the range of 
the eigenvalue distribution does not change over the interface, 

x{q,z) = [ dq'p{q',z). (17) 

J-qO 

Using the chain rule, we have 

dq(x, z) dq(x, z) dx(q, z) 



dz dx dz 

d z f q dq>p(q',z) 



The kinetic term becomes (c.f. [27]) 



(18) 



J p{q,z) 

A. Interface tension at 

For simplicity consider the interface tension at first. At the transition there are both 
order-order interface and order-disorder interfaces. Due to the complicated form of the 
kinetic energy Eq. (19), solving for the full solution of the interface seems to be hopeless. 
However, the presence of the zero modes allows us to show that both interface tensions 
vanish. It is straightforward to construct an interface using the shift mode: 

q = | Q(z + 1 - A/(*)) - 1 + Af{z), 0<x< Af(z); 
q [Q(x-Af(z)) + Af(z), Af(z)<x<l 

with /(— L) = and f(L) = 1 at two boundaries of the interface. A fc-wall interpolating 
two vacua related by Z^- transformation corresponds to A = -j^. 

The idea is to take f(z) — ^ + 1 such that ~ ^- As / is a zero mode, this not change 
the potential energy, while the kinetic energy is supressed by ^. In the limit L — > oo, action 
vanishes as 1/L, and so the interface tension vanishes at T^. 

This does not work for the shift mode. To see that, we need to take a close look at the 
kinetic term. On the interface, the density of eigenvalues is given by 

l + cos27r(g + l-A/(z)), Q(l - Af{z)) - 1 + Af < q < -1/2 + Af(z) 
' ^ l + cos27r(g-A/(z)), -1/2 + Af{z) < q < Q(l - Af{z)) + Af(z). 
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Integrating with respect to q, we obtain 

x = q+ l +±- S in[2vr(g - A/)], Q(l - Af) - 1 + Af < q < Q(l - A/) + Af, (22) 

2 Z7T 

from which it follows 

dq = cos[2n(q-Af)]Af 

dz~~ 1 + cos[27r(g - A/)] ' 1 j 

Plugging the gradient Eq. (23) into Eq. (19), we identify a non-integrable singularity at 
q = 1/2 + A /. Therefore, we conclude that the interface built by the shift mode is ruled out 
by the divergent kinetic energy. 

A second possibility is to build an interface with the 6-mode. The two vacua at the ends 
of the interface are joined through a confining phase in the middle. Defining as the 
inverse function of 

1 b 

x = q + - + —sm(2nq). (24) 

At T = T d , Qi(x) reduces to the distribution for the deconfined phase, and Qo(x) for the 
confined phase. The interface is constructed as 

part I :q = Qb( z )( x ), ~L < z < -H- 

j Qb(z)(x + 1-A)-1 + A, 0<x<A; 
part II :q = < < z < L, (25) 

[ Qb{z){x- A) +A, A<ac<l. 

with the boundary conditions b(±L) = 1, 6(0) = 0. The distributions from part I and part 
II joining at z = are identical and given explicitly by q = x — 1/2. We note that part 
II is a Z(N) transform of part I, flipped in z. We will show below that the kinetic energy 
is not divergent, so that we can apply the previous argument to arrive at the suppression 
in 1/L. Since the path given by the first and second lines of Eq. (25) are related by Zn 
transformations, they necessarily have the same kinetic and potential energies. It is sufficient 
to restrict ourselves to the first line. It is easy to find the gradient 

dq = V sin 2irq 

dz~ 2tt 1 + 6cos27rg' [ ' 

where here and in the following primes denote the derivatives with respect to the argument. 
We see possible singularities at q = ±1/2 from the denominator are cancelled by the nu- 
merator, giving rise to a finite result for the kinetic energy. For a given b, it is not difficult 
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to compute the integral in q by contour integration, 

r iqp(q) m 2 = ( *v 2 - hi ± fc2 < 5 - 3 ^> + 4 (- 1 t^HB . (27 ) 

We checked that the kinetic term is finite in the limit b — > 1. 
B. Interface tension at T > T d 

The order-disorder interface tension is only defined at T d . Above T d) the order-order 
interface tension, equivalent to the 't Hooft loop, is nonzero. However, we were not able to 
compute the order-order interface tension in full generality. 

Let's consider two limiting situations, A — > and T — > T d . We consdier them in turn. 

In the case A — )• 0, the end points of the interface are given by 

}g(x+l-A)-l + A,0<x<A; 
q = Q(x), q = < (28) 
[Q{x-A) + A, A<x<1. 

The eigenvalues between A < x < 1 have change infintesimally, while those between < 
x < A have a finite jump; the latter, however, are suppressed because there are few of them. 
With this in mind, we write down the following path for the interface: 

q = (Q(x + l-A)-l + A-Q(x))g(z) + Q(x), < x < A; (29) 

q = {Q{x-A) + A-Q(x))f{z) + Q{x), A < x < 1. (30) 

The unknown functions g(z) and f(z) interpolate between and 1. They could in principle 
depend on x, which characterizes the change of shape of the eigenvalue density. To leading 
order in A we can ignore this dependence. Now we can work out the potential energy along 
the path 

pA pi pA pi 

V = dxdyV (q(x) - q(y)) + / dxdyV(q(x) - q(y)) + 2 dx dyV(q(y) - q(x)). 

JO J A JO J A 

(31) 

The first term is of order 0(A 3 ) and may be ignored. The second term starts with the 
vacuum potential energy as that at leading order. The third term is of order 0(A). We 
need to know the 0(A) correction of the potential to the vacuum one: 

6V= fdxdy 9V A\Q\x)-Q\y)\f 
Jo d\q{x)-q(y)\ 



+ 2 f dx f dy (V(q(y) + g ) - V(q(y) + ?„ + (1 - 2q )g)) . (32) 
Jo Jo 

Evaluating Eq. (32) for the vacuum solution Eq. (11a), we find the first term vanishes 
identically, while the second term gives 

d 2 



The kinetic energy reads 



5V=^(l-2q )Y(l-g) 2 . (33) 



dzdx [ ^) (34) 



wfth a ?= (2* -IV, 0<.<A; 
^ I A (1 — Q'(x)) f, A<x<l. 



We ignore the contribution from A < x < 1 because it is of order 0(A 2 ). Combining the 
kinetic and potential terms, we have 



J dz(K + 5V) = J dz (A(l - 2q ) V 2 + A^(l - 2q fg 



^ 2 'l~gf). (36) 



We need to minimize the above action with the boundary condition g(—oo) = and g(oo) = 
1. Using a trick (if (x — y) 2 > 0, x 2 + y 2 > 2xy), we can obtain the minimum without solving 
for g: 



j dz(K + 5V)>Aj dz2\j (1 - 2g ) V 2 ^(l - 2q ) V(l - <7) ! 



A/^-^^l^^^^, (37 ) 



which leads to 



« = A^=(l-2g ) 3 . (38) 

Now we look at the other limiting case: T — > We expect the interface to mimic 
the T = Td case, i.e. two vacua joined through a confined distribution. We consider the 
following simple path: 

part I :p = a(z) + b(z)bo cos dq, —qo < q < go, 

-g _ (1 _ A)(l - 2q )f(z) < q < -(1 - 2A)g - (1 - A)(l - 2q )f(z) 



part II :p 



h:i -(1 - 2A)g + A(l - 2q Q )g(z) < g < g + A(l - 2q )g(z) 



2? 



a(z) + b(z) cos d(q + I- A), Q 6 (l - A) - 1 + A < g < g - 1 + A; 
part III :p = ^ , (39) 

o(«) + b(z) cos rf(g + 1 - A), -g + A < g < Q 6 (l - A) + A, 
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We have —L < z < — ^ for part I, — ^ < z < % for part II and ^ < z < L for part III, with 
the following boundary conditions 



b(±L) = 1, b(±L/2) = 0, f(-L/2) = 0, f(L/2) = 1. (40) 

In Eq.(39) g and b are determined by Eqs. (11c) and Eq. (lid), respectively. The 
normalization condition J dqp(q) = 1 forces a(z)go + ^ sin dqo = 1/2, while the tracelessness 
condition J dqp(q)q = sets f(z) = g(z). Furthermore, we have the boundary conditions 
The interface is composed of three parts. Part I makes use of the 6-mode to deform the 
eigenvalue distribution. At T — > T^, b mode is nearly a zero mode, with minimum cost of 
potential energy. Part II is unique to the case T > Tj, with the motion of the eigenvalues 
suppressed by (1 — 2q ). Part II is just a Z N transform of part I, flipped in z. 

Let us look at part I first. The potential energy with respect to the vacuum one has the 
following expansion in | — go (which is an effective expansion in T — Td) 



60 -d 2 ^ ,„ /l 



^ = 1 

90 



il/ = r V L(1 - 6)2 (2-*j + °((2"*j J' <«> 

The kinetic energy can be taken as at T = T^, ignoring higher order corrections Eq. (27). 
With these ingredients, we can already work out the contribution to the interface tension 
from part I. It is worth noting prior doing any computations that 5V ~ (| — go) 2 , K ~ 0(1). 

Next we consider part II. This case is particular easy because of the constant eigenvalue 
density in this part. To leading order, the potential energy and kinetic energy are 

60 -d 2 _60/(l-/)A(l-A)(l2-d 2 A(l-A)) Q - go) , 

T = f 2 A(l-A)Q-g ) 2 . (42) 

We note in part II both potential and kinetic terms are suppressed by (| — go) 2 , leading 
only to contributions of higher order in (T — To) for the interface tension. 

Adding up contributions from all three parts, we obtain a contribution to leading order 
which is twice that of part I. It is not difficult to convince ourselves that a ~ (| — g ) an d 
it is independent of A. 

We will improve the result by considering a more sophisticated ansatz. Note that in Eq. 
(39), we have chosen to turn on the 6-mode and separating /-mode separately. Here we 
consider more general ansatz by turning them on simultaneously. Furthermore, as we learn 
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from the previous example that part II has only subleading contributions, we consider the 
abovementioned modification to part I only. With these in mind, we consider the following 
ansatz 



p = a(z) + b(z) cosd(q + (1 - A)f L (z)), -q - (1 - A)f L < q < Q b (A) - (1 - A)f L 
p = a{z) + b(z) cosd(g - Af R (z)), Q b (A) + Af R <q<q + Af R . (43) 



The tracelessness condition forces f L = f R = f. It is also natural to require / ~ 1 — 2q such 
that the eigenvalues do not overseparate along the path. The introduction of the /-mode 
induce corrections to the potential in Eq. (42). Defining / = / (| — q ) + O ((| — go) 2 ), 
we only need to consider corrections up to the second order in /. It turns out the first 
order in / has a coefficient P(A, — qo) + O ((| — go) 2 ) and the f 2 term has a coefficient 
Q(A, b) + O ((| — qo)). When both are taken into account, we have for the potential 



SV = {5V (b) + P(A,b)f + Q(A,b)f 2 ) | - - u. 



(44) 



0.02 



0.01 - 




FIG. 1. Interface tension at qo = 0.3. Due to the symmetry A -H- 1 — A, we have only shown half 
range of A. The red linear line is for A C 1, while the black dashed line is obtained assuming 
f-2go«l (i.e. T^T d ). 
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where 




SV = 


1 


90 


P = - 


1 




48 



60 - d 2 - 60/(1 - /)A(1 - A) (12 - d 2 A(l - A)) 



(45) 



(6 - 1) (326 + (1 - 4Q b (A) 2 )(-48 + d 2 {l + 8Q b (A) + 12Q 6 (A) 2 - 16AQ 6 (A)))) 



-166(-l + 12Q fc (A) 2 )cosdQ b (A), (46) 
Q = 3846 2 - 486d 2 (1 + 4Q fe (A) 2 - Q 6 (A)(4 - 8 A)) + d 2 48 - 192Q 6 (A) 2 

+ d 2 (-5 - 24Q fe (A) 2 + 48Q 6 (A) 4 + Q fe (A) 3 (32 - 64A) + 16A(1 - A) + 24Q 6 (A)(-1 + 2A)) 
+ 486 (86 + d 2 {l - 4Q 6 (A) 2 )) cos dQ b {A). (47) 

We have also factored out the overall (| — g ) 2 dependence. One important property we 
confirm numerically is that Q > 0. To lower the potential energy, we choose / — As a 



result, 



SV =[6Vo- 



P 2 \ (\ 



go 



(48) 



AQJ \2 

At the same time, the correction also introduces A dependence to the interface tension. 
With the same kinetic energy as at leading order, we obtain the interface tension 



a 



2 J dzvWK~(^-q ). 



(49) 



The final result is obtained numerically and shown in Fig. (1). Both scenarios are included. 

We close this section with the observation that these different scenarios do not have a 
common region in which they are both valid. This suggests the limits T — > Td and A — > 
(or N — > oo) do not commute. Therefore, it is important to evaluate corrections at finite N 
near Td- 



IV. FINITE N CORRECTION NEAR T d 

Corrections at finite N enter where the integrals at N = oo are replaced by discrete sums. 
This can be evaluated with Euler-MacLaurin formula [28] 

(50) 
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where f(i/N) = F(i) and B 2 k are Bernoulli numbers. We wish to shift the argument and the 
lower integration bound from to 0. This reshuffle can be done with a Taylor expansion. 
We end up with the following surprisingly simple expression: 

±F(i) = N f\ m + + f (f ^ {1) _ ^.(oj) . (51) 

i=l " / ° ~ fc=l ^ 

We are not able able to prove Eq. (51), which requires a recursion relation among Bernoulli 
numbers. We have verified explicitly that it holds for the first few terms in the expansion. 
Applying Eq. (51) to the sums, we have 

N 



N = J2 1 = N 

Jo 



dx (52) 

= f q = N f 1 dxq + S^M + ... (53) 
Jo 2 

V=± V M = N> f d.dyV( q ( X ) - q (y)) + 2N f Z gW) Z V JM Z «(")) 

+ -.., (54) 

where • • • denote terms higher order in 4> expansion. As before, we define the eigenvalue 



density p(q) = ^f, in terms of which Eq. (52) becomes 



dq 



1+ 



1 = / dqp(q) (55) 

Jq- 

= P dqp(q)q + + • • • (56) 

V = N 2 [" dqdq'p(q)p(q')V(q - q') + N f + dqp(q) (V(q + - q) - V(q - <?_)) + • • • , (57) 

with g_ = g(0) and g + = q(l). For the infinite N eigenvalue distribution, p(q) = p(—q) and 
q(0) + g(l) = 0. The correction at next to leading order to V vanishes identically. However 
this need not be true for the corrected distribution in ^. 

To determine the new eigenvalue distribution, we vary V with respect to p, subject to 
the usual normalization and tracelessness constraints in Eq. (55). This gives 

2 / dq'p{q')V{q - g') + - (V(q + - q) - V(q - </_)) + A a + X 2 q = (58a) 
dqp{q) = 1 (58b) 
^ dqp(q)q+?±^ = 0. (58c) 
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We follow the method of [8] in solving for p. Taking derivative with respect to q four times, 
we arrive at 



p"(q) + d\p(q)-l) = 0. (59) 

The general solution is given by 

p(q) = 1 + b cos dq + csmdq. (60) 

Plugging it into Eq. (58a), we find the result is organized as a fourth order polynomial 
in q. The coefficients of q and the constant term can always be set to zero by choice of 
the Lagrange multipliers Ai and A2. The remain coefficients of g 4 , g 3 and q 2 give three 
independent equations. Together with Eq. (58b) and Eq. (58c), we have in total five 
equations, to be satisfied by four constants b, c and q±. It turns out when four of the 
equations are satisfied, the fifth automatically holds. In practice, we solve for cos dq± and 
smdq± in terms of the constants b, c and q±. Defining q s = q_ + q + and qd = q + — g_, we 
can first find an equation for q^ 

2 6d(l-q d ) K ' 

It is easy to see qa = 2qo as defined in Eq. (11c) is free of corrections in 1/N, and that 

l_ -12-2rf 2 g (l-2g ) 2 
qs ~ N rf 2 (l-2g ) 2 1 ] 

(144 + Yld\\ - 2g ) 2 + d\\ - 2g ) 4 ) (36 + d 2 N\l - 2g ) 4 ) _ 

d?{l-2q y 1 + j 1 } 

The first shows that q s is suppressed in 1/iV. The second combined with 

tan^= 1 - C °;^ + ^ ) , (64) 
2 smd(q.+q + ) ' v ; 

can be used to determine b and c. To have a consistent i expansion, we need to have 

^ = ^o + ^ + ---,c=| + .... (65) 

At leading order 60 agrees with Eq. (lid). c\ is 

ci = -6 d go- (66) 

We have obtained Eqs. (60), (61), (62), (65) and (66) as the new eigenvalue distribution 
up to order j=. The appearance of sm(dq) in Eq. (60) looks worrisome at first sight: 
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the Polyakov loop can have a nonvanishing imaginary part, violating charge conjugation 
symmetry. The cause of the strange behavior lies in our expansion. We have taken Eq. 
(52) literally as the 4 expanison and used it to solve for the eigenvalue distribution p to the 
first nontrivial order. However, the result contains constants which organize themselves in 

expanison, making our naive expansion inconsistent. A consistent expansion has to take 
into account the h factors in the constants. 

Using the experience from before, we expect that to order i, the constants have the 
following form: 



7 7 C l O QSI 

c = iV' qd= g °' Qs = iV 



(67) 



Plugging Eq. (67) into Eqs. (58a), (58b), (58c) and keeping terms up to order we find, 
surprisingly, that the resulting equations can be solved by 

Ci = ' qsi = - d *(i U 2 qo r (68) 

The term sin(dg) naturally drops out, and we find explicitly that the imaginary part of all 
Polyakov loops vanish to order 4. We summarize the eigenvalue distribution as follows 



Thus we have 



p = 1 + b cos dq, q- < q < q + 
12 

q - +q+ = -d*N(i-2q r q+ ~ q - 



-(q + -q )d 2 N(l-2q ) 2 /6 = l. 



2q . 



(69) 
(70) 

(71) 



In Fig. (2), we compare this combination with the numerical results for the maximal eigen- 
value, g max , which at high N is approximated by q + , for different number of colors and 




1.2 1.4 1.6 1.8 2.0 

d/(2jt) 



FIG. 2. The maximal eigenvalue, q max , as a function of d for different number of colors, N. 
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different d, that is temperature. Nice agreement between analytic expression Eq. (71) and 
numerical results is seen. From Eq. (69) and our numerical simulations we also see that 
when N(l — 2q ) 2 ~ 1, the large N expansion breaks down. We numerically computed that, 
at d = 2tt, i.e. T = T^, we have q ma ,x — Qo ~ l/ViV, in contrast to the behaviour of the 
maximal eigenvalue in the confined and the deconfined phase q ma x — Qo ~ 1/N. 

With the correct eigenvalue distribution, we can proceed to evaluate the correction to 
the potential. We find the correction vanishes to order j^. To obtain the correction to the 
next order, we need to find p at the corresponding order. We will not do it at this time, but 
use Eq. (69) to give an estimate of the correction. We need the following formula for the 
potential 

rg+ 

N 



V 



q+ ] 
dqdq'p(q)p(q')V(q - q) + 

q_ - ' Jq 



dqp(q) (V{q + - q) - V(q - q_)) 



V{q+-q-) B 2 
2N 2 N 2 



(72) 



Taking into account 4 correction to q s , we have i contribution from all three naive orders. 
We will not spell out the detail of the calculation, but only list the final results: 

,3\ 



t/ _ 180 - 5rf 4 g0 2 (l - 2g ) 5 + 24d 2 q 2 (-5 + 20q - 30g 2 + 16^, 
Vi — 



v 2 



30d 2 N 2 (l - 2q f 
g (l-2g )(6- Ao(l-2g )) 
3A^ 2 

go(l-2go)(6-Ao(l-2g )) 
6iV 2 



(73) 

(74) 
(75) 



where Vi, V 2 and V 3 denote contributions from naive order 1, 4> and respectively. We 



note the large N expansion breaks down as iV 2 (l — 2q ) 



, N 

1. The bulk thermodynamic 



quantity has a less stringent criterion than the eigenvalue distribution. 




-0.4 -0.2 o 



FIG. 3. The non-equilibrium potential U(L) = V(L) - V(L = 0). 
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FIG. 4. The maximum of the potential U max for different N. 

We were able to analytically estimate corrections to the potential at the minumum only. 
Numerically it is possible to go beyond this and compute the potential as a function of the 
Polyakov loop L. In Fig. (3), the potential normalized by N 2 — 1 is shown as a function of 
L for different number of colors. It is remarkable, that the potential flattens from the side 
of the confinement point = 0. The potential maximum characterizes the order-disorder 
and order-order interface tensions and it falls like 1/N 2 for large N, which is in qualitatively 
agreement with large N expansion, as demonstrated in Fig. (4). 
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V. CONCLUSIONS 

The matrix model studied here is clearly only one of many possible matrix models. Its 
advantage is that it can be solved analytically at infinite N, and numerically at finite N. 
In this model all interface tensions vanish at the deconfining phase transition. Since the 
transition is of first order, this would be striking evidence that it is an unusual transition, 
perhaps of the Gross- Witten-Wadia type. 

Presently, numerical simulations of SU(N) gauge theories can only be carried out at 
relatively small N, N < 10. For two colors the order-disorder interface tension vanishes, as 
the transition is of second order. For three colors one the order-disorder interface tension is 
relatively small [25]. This presumably reflects that the transition for three colors is weakly 
first order, because of its proximity to the second order transition for two colors. This leads 
one to expect that as iV increases, that the order-disorder interface tension, divided by N 2 , 
increases monotonically from N = 3, and becomes constant at infinite N. 

The numerical solution of the matrix model indicates the contrary, that the order-disorder 
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interface tension, divided by iV 2 , behaves non-monotically with N. From Fig. (4), the 
barrier of the potential, suitably normalized, increases from N = 3 to N = 5, and then 
slowly decreases as N increases further. 

Such non-monotonic behavior in N is unexpected, and could well just be an artifact of 
the model. This could be settled by numerical studies on the lattice of the order-disorder 
interface tension for moderate values of iV ~ 5. It might even provide hints of a Gross- 
Witten-Wadia transition at infinite N. 
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